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Abstract 

In this paper, we review the so-called Topfer algorithm that allows us to 
find a non-iterative numerical solution of the Blasius problem, by solving a 
related initial value problem and applying a scaling transformation. More- 
over, we remark that the applicability of this algorithm can be extended to 
any given problem, provided that the governing equation and the initial con- 
ditions are invariant under a scaling group of point transformations and that 
the asymptotic boundary condition is non-homogeneous. Then, we describe 
an iterative extension of Topfer's algorithm that can be applied to a general 
class of problems. Finally, we solve the Falkner-Skan model, for values of 
the parameter where multiple solutions are admitted, and report original nu- 
merical results, in particular data related to the famous reverse flow solutions 
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by Stewartson. The numerical data obtained by the extended algorithm are 
in good agreement with those obtained in previous studies. 
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1 Introduction 

At the beginning of the last century, Prandtl [|43l put the foundations of boundary- 
layer theory providing the basis for the unification of two sciences, which at that 
time seemed incompatible: namely, theoretical hydrodynamics and hydraulics. 
The main application of boundary-layer theory are devoted to the calculation of 
the skin-friction drag acting on a body moving through a fluid, for example the 
drag of: an airplane wing, a turbine blade, or a complete ship (see Schlichting 
and Gersten ATI ). In this new century, due to the increasing number of applica- 
tions of microelectronics devices, boundary-layer theory has found a renewal of 
interest within the study of gas and liquid flows at the micro- scale regime (see 
for instance Gad el Hak [[HI or Martin and Boyd [40]). Recently, Boyd [[lOl has 
used the problem considered by Prandtl as an example where some good analysis 
allowed researchers of the past to solve problems, governed by partial differential 
equations, that might be otherwise impossible to face before the computer inven- 
tion. In this context, Prandtl, Blasius and Topfer work is still modem now as it 
was more than one hundred years ago. 

For an exhaustive derivation of the complete boundary layer equations govem- 
ing general compressible flows, the interested reader is referred to Stewartson [|50l 
or Schlichting and Gersten [|47l . Within boundary-layer theory, Blasius problem 



2 



dHIl is given by: 

where t] and / are suitable similarity variables. This is a two point boundary value 
problem (BVP) obtained for the model describing the steady plane fluid flow past 
a thin flat plate. 

Boyd has pointed out how this particular problem of boundary -layer theory 
has arose interest of prominent scientists, like H. Weyl, J. von Neumann, M. Van 
Dyke, etc.; see Table 1 in |i9|. The main reason for this interest is due to the hope 
that any approach developed for this epitome can be extended to more difficult 
hydrodynamics problems. 

Blasius main interest was to compute, without worrying about existence or 
uniqueness of its BVP solution, the value of 



^ = 0(0), (1.2) 



i. e., the skin-friction coefficient. In order to compute this value, Blasius used a 
formal series solution around 7] = and an asymptotic expansions for large values 
of 7], adjusting the constant A so as to connect both expansions in an intermediate 
region. In this way, Blasius obtained the (wrong) bounds 0.3315 < A < 0.33175. 

In 1912, Topfer [l52l revised the work by Blasius and solved numerically the 
Blasius equation with initial conditions 



Using hand computations with the classical fourth order Runge-Kutta method and 
a suitable scaling invariance, he arrived, without detailing his computations, at the 
value A ^ 0.33206, contradicting the bounds obtained by Blasius. 



Thereafter, the quest for a good approximation of A became a main concern. 
This is seldom the case for the most important problems of applied mathemat- 
ics: the first study is usually devoted to find a method to solve a given problem, 
and once the problem is solved the attention is turned to know how accurate is 
the computed solution and whether there are different methods that can provide 
a solution with less effort or more accuracy. Using a power series, Bairstow [jVl 
reports A ~ 0.335, whereas Goldstein [|35l obtains A ~ 0.332; moreover, using 
a finite difference method, Falkner fl9] computes A ~ 0.3325765, and Howarth 
[l37ll yields A ~ 0.332057. Fazio If24l . using a free boundary formulation of the 
Blasius problem, computes A ~ 0.332057336215. Boyd flUl uses Topfer's algo- 
rithm to obtain the accurate value A ~ 0.33205733621519630. By the Adomain's 
decomposition method, Abbasbandy [TT| finds A ~ 0.333329, whereas a varia- 
tional iteration method with Pade approximants allows Wazwaz [|55l to calculate 
the value A ~ 0.3732905625. Tajvidi et al. [51J apply modified rational Legendre 
functions to get a value of A 0.33209. 

Also the Crocco formulation [fTTl can be applied in order to compute the value 
of skin-friction coefficient, we can apply . For instance, Vajravelu et al. [|53l 
use the Runge-Kutta method and a shooting technique to solve numerically the 
Crocco formulation and obtain the value A ~ 0.3322, whereas Callegari and Fried- 
man [fT2l reformulate the Blasius problem in terms of the Crocco variables, show 
that this problem has an analytical solution, and compute the following bounds: 
0.332055 < A < 0.33207. 

Within scaling invariance theory, the transformation of a BVP into an initial 
value problem (IVP) due to Topfer is a consequence of the invariance of Blasius 
equation and the two boundary conditions at 7] = (at the plate) with respect to 
the scaling group of transformations 

r=A-«/, T?*=A«7], (1.4) 
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where a is a non-zero parameter. Using this transformation, a simple existence 
and uniqueness theorem was given by J. Serrin, as reported by Meyer [|4T1 pp. 
104-105]. This scaling invariance is essential also to the error analysis of the 
truncated boundary formulation for Blasius problem due to Rubel 1441 . A more 
complex proof of the existence and uniqueness of the Blasius problem solution 
was given before by Weyl [|56l. who also proved that the solution has a positive 
second order derivative that is monotone decreasing on [0, °°) and approaches to 
zero as rj goes to infinity. 

Our main goal here is to show how to solve numerically Blasius problem and 
similar problems in boundary layer theory, by initial value methods defined within 
scaling invariance theory. These methods are referred to as numerical transforma- 
tion methods (TMs). As pointed out by Na Il42l Chapters 7-9], usually a given, 
even simple, extension of the Blasius problem cannot be solved by Topfer algo- 
rithm. Therefore, in order to extend the applicability of this non-iteratve transfor- 
mation method (ITM) an iterative version has been developed in [|26ll27]| . 

In particular, we apply the ITM to the Falkner-Skan equation with relevant 
boundary conditions 



drj^ drj 



drj 





(1.5) 



where / and r] are appropriate similarity variables and /3 is a parameter. This 
set is called the Falkner-Skan model, after the names of two English mathemati- 
cians who first studied it [|20l. As pointed out by Na [l42l pp. 146-147], if /3 7^ 0, 
the BVP (11.51) cannot be solved by a non-ITM. Indeed, the goveming differen- 



tial equation in (11.51) is not invariant with respect to any scaling group of point 
transformations. 

The existence and uniqueness question for the problem (|1.51 ) is really a com- 



plex matter. Assuming that /3 > and under the restriction < ^ < 1> known as 
normal flow condition, Hartree [l36l| and Stewartson [|49l proved that the problem 
(|1.5I) has a unique solution, whose first derivative tends to one exponentially. Cop- 
pel lfT4l and Craven and Peletier [[T5ll pointed out that the above restriction on the 
first derivative can be omitted when < /3 < 1 . Weyl proved, in [|56]| . that for each 
value of the parameter j8 there exists a physical solution with positive monotone 
decreasing, in [0,oo), second derivative that approaches zero as the independent 
variables goes to infinity. In the case /3 > 1, the Falkner-Skan model loses the 
uniqueness and a hierarchy of solutions with reverse flow exists. In fact, for /3 > 1 
Craven and Peletier [[T6l computed solutions for which ^ < for some value of 
7] . In each of these solutions the velocity approaches its limit exponentially in rj . 
As mentioned before, the term normal flow indicates that the flow velocity has a 
unique direction, and instead, reverse flow means that the velocity is both positive 
and negative in the integration interval. 

The considered problem has also multiple solutions for /3inin < /3 < 0, as re- 
ported by Veldman and van de Vooren [|54l . with the minimum value of /3 given 
by 

/3min = -0.1988... . (1.6) 

In this range there exist two physical solutions, one for normal flow and one for 
reverse flow. For /3 = jSmin only one solution exists. Finally, for /3 < /3niin the 
problem has no solution at all. Our interest here is to apply the ITM to the range 
of /3 where multiple solutions are admitted, in particular, our interest is to get 
numerically the famous solutions of Stewartson [|49ll50l . The obtained results are 
original and, as we shall see in the following sections, are in agreement with those 
available in literature. 

The first computational treatment of the Falkner-Skan model is due to Hartree 
[|36ll . Cebeci and Keller 111311 apply shooting and parallel shooting methods requir- 
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ing asymptotic boundary condition to be imposed at a changing unknown bound- 
ary in the computation process. As a result, they report convergence difficulties, 
which can be avoided by moving towards more complicated methods. Moreover, 
to guarantee reasonable accuracy, they are forced to use a small enough step-size 
and extensive computation for the solution of the IVPs. Na [42, pp. 280-286] 
describes the application of invariant imbedding. A modified shooting method [0 
and finite-difference methods p3l,l4l| for this problem are presented by Asaithambi. 
Kuo If39ll uses a differential transformation method, which obtains a series solu- 
tion of the Falkner-Skan equation. Sher and Yakhot [|481 define a new approach 
to solve this problem by shooting from infinity, using some simple analysis of the 
asymptotic behaviour of the solution at infinity. Asaithambi O proposes a faster 
shooting method by using recursive evaluation of Taylor coefficients. Zhang and 
Chen [57] investigate a modification of the shooting method, where the computa- 
tion of the Jacobian matrix is obtained by solving two IVPs. A Galerkin-Laguerre 
spectral method is defined by Auteri and Quartapelle [0. 



2 Topfer transformation 



In order to clarify Topfer [|521 derivation of a transformation of variables that re- 
duces the BVP (11.11) into an IVP we will consider the derivation of the series 
expansion solution. Of course, some of the coefficients of the series can be eval- 
uated by imposing the boundary conditions at r] = 0. Moreover, for the missing 
initial condition, we set 

where A is a non-zero constant. So that, we look for a series solution defined as 

oo 

/(r?) = y^'+lQT]" 

n=3 
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where the coefficients A and C„, for n = 3, 4, . . . , are constants to be determined. 

In fact, the boundary values at the plate surface, at r] = 0, require that Cq = Ci = 0, 
and we also have C2 = A/2 by the definition of A. Now, we substitute this series 
expansion into the governing differential equation, whereupon we find 



£n(n-l)(n-2)C„T]"-3+ 

n=3 



A + J^n(n-1)C„T] 

n=3 



n-2 



= 



or in expanded form 

[3-2-C3] + [4.3-2-C4]t] + 

+ 



5.4.3-C5+i-24T 

2 2 2 



6-5-4-C6+i-2-^C3 + ^-^-3-2-C3 

° 2 2 2 2 



T]^ + --- = 0. 



According to a standard approach we have to require that all coefficients of the 
powers of J] to be zero. It is a simple matter to compute the coefficients of the 
series expansion in terms of A : 

A2 



€3 = €4 = 0, 

C6=C7 = 0, 

C9 — Cio — , 
and so on ... 



C5 



C8 = ll 



2-5! 
A3 

22-8! 



Cii = -375 



23-11! 



The solution can be written as 

At]2_AV ll-A^T]^ 375 -A^T?" 
2-5! ^ 22-8! 23.11! 

where the only unknown constant is A. In principle, A can be determined by 
imposing the boundary condition at the second point, but in this case this cannot 
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be done because the left boundary condition is given at infinity. However, by 

modifying the powers of X we can rewrite the series expansion as 

[X'l^n^ (^1/377)^ 11-(ai/377)^ 375-(Ai/3r7y' 
^ 2 2-5! ^ 22-8! 23.11! 

which suggests a transformation of the form 

f = X-^l^f^ r]* = X^I^'{]. (2.1) 

In the new variables the series expansion becomes 

ri*2 ri*5 u.ri*s 375.17*11 
~^~2^^ 22-8! 23-11! 

which does not depend on X. We notice that the governing differential equation 
and the initial conditions at the free surface, at r] = 0, are left invariant by the 
new variables defined above. Moreover, the first and second order derivatives 
transform in the following way 

d-q* d-q ' dq*^ dr]^ ' 

As a consequence of the definition of A we have 

and this explain why in these variables the series expansion does not depend on A. 
Furthermore, the value of X can be found provided that we have an approximation 

df\ 



for ^ (00). In fact, by the above relation we get 



X = 



df* 



(00) 



-3/2 



(2.2) 



From a numerical viewpoint, EVPs must be solved within the computational 
domain simultaneously (a "stationary" problem), whereas IVPs can be solved by 
a stepwise procedure (an "evolution" problem). Somehow, numerically, IVPs are 
simpler than EVPs. 



2.1 Topfer algorithm 

Let us list the steps necessary to solve the Blasius problem by the Topfer algo- 
rithm: 



1 . we solve the IVP 



/*(0) = ^(0)=0, ^^{0) = l 



(2.3) 



and, in particular, get an approximation for ^^(°°) in order to compute A 



by equation (12.21) : 



df d^f 

2. we obtain /( 7]), "i— (t?), and -r^(^) by the inverse transformation of (12. It . 



In this way we define a non-ITM. 

Indeed, Topfer solved the IVP for the Blasius equation once. At large but finite 
T]*, ordered so that r\* < we can compute by equation (12.21) the correspond- 
ing Xj. If two subsequent values of Xj agree within a specified accuracy, then A 
is approximately equal to the common value of the Xj, otherwise, we can march 
to a larger value of X] and try again. Using the classical fourth order Runge-Kutta 
method (see Butcher [[TTl p. 166]) and a grid step A?]* = 0.1 Topfer was able to 
determine X with an error less than 10^^. He used the two truncated boundaries 
7)1 = A and r/l = 6. We reproduce Topfer computations in figure H] In fact, in 
figure [T] we plot the numerical solutions obtained by Topfer' s algorithm defined 
above. We notice that the top and the bottom frames of this figure show the solu- 
tions of the IVP (O and of the BVP (fTTl) . 
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3 Extension of Topfer algorithm 



The applicability of a non-lTM to the Blasius problem is a consequence of its 
partial invariance with respect to the transformation (12.11) : the asymptotic bound- 
ary condition is not invariant. The non-iterative algorithm can be extended, to a 
given problem in boundary layer theory, provided that the goveming equation and 
the initial conditions are invariant under a scaling group of point transformations 
and the asymptotic boundary condition is non-homogeneous. Several problems 
in boundary-layer theory lack this kind of invariance and cannot be solved by 
non-ITMs [42, Chapters 7-9]. To overcome this drawback, we can introduce an 
iterative extension of the algorithm. The main idea for the new algorithm is to 
modify the original problem, by introducing a numerical parameter h, and to re- 
quire the invariance of the modified problem with respect to an extended scaling 
group involving h; see |l26l|27|| for details. 

3.1 The iterative transformation method 

In order to define the ITM, let us consider the class of EVPs defined by 



Jt]3 ^y^^^^dii'dv^ ^^^^ 

m=a ^(0)=^, ^(^)^^ 

where a, b and c 7^ are given constants. We modify the class of problems (13.11 ) 
by introducing a numerical parameter h as follows 

Jr]3 ^^y' ^''^ ^'^ dy]'"" dr]^) ^^^^ 

f(0)=h^/''a ■^(0)=/z(i-^)/^Z>, ^(r])^c as t]^oo. 
dr] drj 

It is worth noticing that the extended problem (13.21) reduces to the original prob- 
lem (|3.1I) for h — 1. Moreover, the extended problem (|3.21 ) is partially invariant, 
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the asymptotic boundary condition is not invariant, with respect to the extended 
scaling group of transformations 

r=A/, T]* = A^T], /i*=A^/7, (3.3) 

with 8^1 and a 7^ 0. Therefore, to find a solution of the given BVP means to 
find a zero of the so-called transformation function 

r(r) = A^'^r-l , (3.4) 

where the group parameter A is defined by the formula 

1/(1-5) 

(3.5) 

and to this end we can use a root-finder method. Let us notice that A and the 
transformation function are defined implicitly by the solution of the IVP 

_ , *(l-35)/cT f,*-iS/cy) .-l/a A5-l)lad£_ A28-l)/a (ff* \ 

2 * ^^-^^ 

r (0) = h)'/% ^(0) = h*^'-^^/% , ±11(0) =d, 

where J is a parameter fixed by the user. In particular, we are interested to com- 
pute ^^{t\V)^ an approximation of the asymptotic value ^^{°°), which is used in 
the definition of A . 

We set the values of d and a and follow the steps: 

1 . we apply a root-finder method to compute a sequence h* , for j = 0, 1 , 2, . . . , 
. Two sequences Aj and T{h*) for j = 0, 1, 2, ... , are defined by equation 
(|3.5I) and (13.41) . respectively. 

2. a suitable convergence criterion should be used to verify whether r(/z*) — )• 
as j — )■ 00. In this way we find the correct value of h* that transforms into 
h=\.lf this is the case, then A; converges to the correct value of A in the 
same limit. 
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df^ 
dri* 



3. once the correct value of A has been found, the solution of the original 
problem can be obtained by rescaling. In particular, we have that 

^(0)=A2^-i^(0). 

In this way we define an ITM. 

In the next sub-section we apply the above iterative extension of Topfer algo- 
rithm to the Falkner-Skan model. 



3.2 The Falkner-Skan model 



In order to apply an ITM to (11.51) we have to embed it to a modified model and 
require the invariance of this last model with respect to an extended scaling group 
of transformations. This can be done in several ways that are all equivalent. In 
fact, the modified model can be written as 



dr\^ dr\ 



drj 







(3.7) 



and the related extended scaling group is given by 

r = A/, T]*=A-ir], r=A^/i, (3.8) 

where a is a parameter. In the following we set a = 4; for the choice a = 8 
see [|26l . In [|26B, a free boundary formulation of the Falkner-Skan model was 
considered and numerical results were computed for the Homann flow (/3 = 1/2) 
as well as for the Hiemenz flow (/3 = 1). 

From a numerical point of view the request to evaluate ^(°°) cannot be ful- 
filled. Several strategies have been proposed in order to provide an approximation 
of this value. The simplest and widely used one is to introduce, instead of infinity, 
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a suitable truncated boundary. A recent successful way to deal with such a issue 
is to reformulate the considered problem as a free B VP [|24] l26l l27l ; for a survey 
on this topic see ||3TI . Recently, Zhang and Chen [[57ll have used a free bound- 
ary formulation to compute the normal flow solutions of the Falkner-Skan model 
in the full range /3inin < /3 < 40. They applied a modified Newton's method to 
compute both the initial velocity and the free boundary. For the sake of simplic- 
ity, we do not use the free boundary approach but, following Topfer, we use some 
preliminary computational tests to find a suitable value for the truncated boundary. 
At each iteration of the ITM we have to solve the IVP 



£_ -L /■* ■>_ I ft 



dri* 





(3.9) 



/*(0) = 4^(0)=0, ^(0) = ±1. 

Tables [U and [21 list the numerical iterations obtained for a sample value of /3. We 
notice that we solve an IVP governed by a different differential equation for each 
iteration because the Falkner-Skan equation is not invariant under every scaling 
group of point transformation. We have chosen /3 = —0.01 since, in this case, 
the missing initial conditions for the normal and reverse flows are not symmetric 
with respect to the /3 axis. The data listed in tables [T] and [2] have been obtained by 
solving the modified Falkner-Skan model onrj* E [0, 20] by setting 

d^ f* 

respectively. In both cases, we achieved convergence of the numerical results 
within seven iterations. Let us now investigate the behaviour of the transformation 
function. Figure |2] shows r(/i*) with respect to h* for the two cases reported in 
these tables. The unique zero of the transformation function is marked by a circle. 
It is worth noticing that the same scale has been used for both axes. As it is easily 
seen, in both cases, we have a monotone increasing function. We notice on the 
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left frame, corresponding to a normal flow, that the tangent to the F function at 
its unique zero and the h* axis define a large angle. This is important from a 
numerical viewpoint because in such a case we face a well-conditioned problem. 
On the other hand, this is not the case for the function plotted on the right frame of 
the same figure. The meaning is clear, reverse flow solutions are more challenging 
to compute than normal flow ones. Therefore, one has to put some care when 
choosing the convergence criteria for the root-finder method. 

Figure |3] shows the results of the two numerical solutions for a different value 
of /3, namely /3 = —0.15. In the top frame we have the normal flow and in the 
bottom frame we display the reverse flow solution. In both cases the solutions 
were computed by introducing a truncated boundary and solving the IVP in the 
starred variables on rj* E [0,20] with = \, h\ = 5 in the top frame and = 
15, h\ = 25 in the bottom frame. In this case, we achieved convergence of the 
numerical results within eight and seven iterations, respectively. For the sake of 
clarity, we omit to plot the solutions in the starred variables computed during the 
iterations. Moreover, we display only 7] e [0, 10]. 

In table [3] we report data concerning the computational cost of the IVP solver 
for the case /3 = —0. 15. Here steps is the number of successful steps, failed is the 
number of failed steps, and evaluations are the calculated function evaluations. 

As far as the reverse flow solutions are concerned, in table 5] we compare the 
missing initial condition computed by the ITM for several values of j8 with results 
available in literature. The agreement is really good. It is remarkable that among 
the studies quoted in the introduction only a few report data related to the reverse 
flow solutions. 

In figure [5] we plot the behaviour of missed initial condition versus j8 . The 
solution found by the data in table [2] is plotted in this figure, but not the one found 
in table [Hbecause this is very close to the Blasius solution. A good initial choice 
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of the initial iterates of h*, for a given value of j8, is obtained by employing values 
close to the one used in a successful attempt made for a close value of /3. It 
is interesting to note that, for values of /3 < /3min the ITM continued to iterate 
endlessly, whatever set of starting values for h* are selected. 

Our extended algorithm has shown a kind of robustness because it is able to 
get convergence even when, for a chosen value of h*, the IVP solver stops before 
arriving at the selected truncated boundary getting a wrong value of r{h*) = — 1. 
On the other hand, the secant method gives an overflow error when this happens 
for two successive iterate of h*. 

The value of /3niin, corresponding to a separation point at T] = 0, can be found 
by the ITM by considering /3 as a continuation parameter. As we have seen, for 
J^min < /3 < two solutions are available: a positive and a negative skin-friction 
coefficient, the missing initial condition, providing a normal and reverse flow so- 
lution. For instance, when /3 = —0. 1988 we get for the missing initial conditions 
the values 0.005221 and —0.005158, respectively. Starting from this value of /3 
we can reduce it gradually and check whether the two missing initial conditions, 

,2 r 

the positive and negative values of ^^(0), converge to zero. Soon, we realize 
that we are forced the use the ITM to its natural limit. In fact, we are trying to 
get a skin-friction coefficient close to zero rescaling a fixed non-zero value, plus 
or minus one in our case. Anyway, when /3 = —0.198837723795 we found the 
skin-friction coefficients 6.61D — 06 and — 6.61D — 06 with 20 and 24 iterations, 
respectively. Finally, we have noticed that, as far the guest for this limiting value 
of /3 is concerned, we are allowed to reduce the chosen truncated boundary value, 
and for /3 = —0.198837723795 this truncated boundary was set equal to one, i.e. 
all IVP was solved on [0, 1]. 

In figure |4] we plot the unique solution for the limiting value /3min, where ji^nin 
is given by equation (11.61) . As it easily seen this is a normal flow solution. 
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The results reported so far have been found by a variable order adaptive multi- 
step IVP solver that was coupled up the simple secant method. The adaptive 
solver uses a relative and an absolute error tolerance, for each component of the 
numerical solution, both equal to then to the minus six. As well known, the secant 
method is convergent provided that two initial iterates sufficiently close to the root 
are used, and its convergence is super-linear with an order of convergence equal to 
( 1 -t- v^)/2. As far as a termination criterion for the secant method is concerned, 
we enforced the conditions 

|r(/ip|<Tol and < TolR|/i*| +T0IA , (3.10) 

with Tol = TolR = TolA = ID — 06. All computations were carried out on a 1.79 
GHz AMD Turion processor with 1 GB of RAM, and for each case the execution 
time was few seconds. 

4 Conclusions 

In this paper we have shown how the original treatment of the Blasius problem 
due to Topfer can be extended to more complex problems of boundary layer the- 
ory. Our main concern was to solve numerically the Blasius problem, and similar 
problems in boundary layer theory, by initial value methods derived within scaling 
invariance theory. In particular, we consider the Falkner-Skan model, for values of 
the parameter where multiple solutions are admitted, and report and compare the 
obtained numerical results, in particular data related to the famous reverse flow 
solutions of Stewartson. 

As mentioned in the introduction, the Falkner-Skan model is a classical ex- 
ample where the simple shooting method cannot be applied satisfactorily for all 
values of the parameter j8. In fact, by applying the usual shooting method to the 
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Falkner-Skan model, one frequently obtains floating-point overflows in the calcu- 
lations (see, for instance, Asaithambi [|2l|5l). A further difficulty is that the initial 
estimate of the missing initial condition must occasionally be very close to the 
exact value in order to get convergence, cf. Cebeci and Keller [1131 

The ITM has the same conceptual simplicity of the simple shooting method. 
It is an initial value method even if we solve a different model for each itera- 
tion when the goveming differential equation is not invariant under every scaling 
group of point transformation. Its versatility has been shown by solving several 
problems of interest: free boundary problems [|33l l23l ^ ^9|, a hyperbolic mov- 
ing boundary problem [|25l . the Homann and the Hiemenz flows govemed by 
the Falkner-Skan equation in [26|, one-dimensional parabolic moving boundary 
problems [|30l . two variants of the Blasius problem [321 . namely: a boundary 
layer problem over moving surfaces, studied first by Klemp and Acrivos [|38l . 
and a boundary layer problem with slip boundary condition, that has found ap- 
plication to the study of gas and liquid flows at the micro-scale regime [fT8l l40l . 
parabolic problems on unbounded domains [|34l and, recently, see the preprints: 
[|22l[ parabolic moving boundary problems, and [1271 an interesting problem in 
boundary layer theory: the so-called Sakiadis problem [|45ll46l . 

In particular, in BOl[ the ITM is used to solve the sequence of free bound- 
ary problems obtained by a semi-discretization of ID parabolic moving boundary 
problems. In [|22l a class of parabolic moving boundary problems is transformed 
to free boundary problems govemed by ordinary differential equations that can 
be solved by the ITM. And in [f34l[ a free boundary formulation for the reduced 
similarity models, that can be solved by the ITM, is used in order to propose a 
moving boundary formulation for parabolic problems on unbounded domains. 
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Figure 1: Blasius solution solved with Topfer's algorithm. Top frame: for 7]* G 
[0,4] we get Ai = 0.3329124105. Bottom frame: for 7]* e [0,6] we find A2 = 
0.3320575595. 
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h* h* 

Figure 2: Two cases of the r{h*) function: left and right frames are related to 
normal and reverse flow solutions, respectively. 



d f* 

Table 1: Iterations for j8 = —0.01 with - — ^(0) = 1. Here and in the following 



drj 



the D — k= 10 means a double precision arithmetic. 



1^1 



(0) 



5. 



0.631459 



0.431723 



1 10. 



1.791425 



0.384034 



2 2.278111 -0.182888 3.389602 0.454658 



2.993420 0.0465208 0.238960 0.454658 



4 2.848366 9.5D-04 0.050925 0.456418 



2.845340 -5.0D-06 0.001064 0.456455 



2.845356 6. ID -08 5.6D-06 0.456455 



2.845355 7.3D-10 6.7D-08 0.456455 
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Figure 3: Normal and reverse flow solutions to Falkner-Skan model for /3 
The symbols • denote values of /(r/). 
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Figure 4: Numerical solutions to Falkner-Skan model for /3 = —0.1988376. We 
(fif 

notice that ^-^(0) = and values of f{T\) are marked by •. 













Blasius 
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Figure 5: Missing initial conditions to Falkner-Skan model for several values of /3 . 
Positive values determine normal flow, and instead, negative values define reverse 
flow solutions. 
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(f- f* 

Table 2: Iterations for B = -0.01 with ^^(0) = -1. 



Y{h)) 



(0) 



75. 



0.731890 



-0.059237 



1 150. 



5.263092 



-0.092368 



2 62.885833 -0.443040 1.385275 -0.028870 



3 69.649620 0.181067 0.097112 -0.046991 



4 67.687299 -0.011297 0.028991 



-0.042016 



5 67.802542 -2. ID -04 0.001700 



-0.042315 



6 67.804749 2.8D-07 3.3D-05 -0.042321 



7 67.804746 7.9D-10 4.3D-08 -0.042321 
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Table 3: Efficiency of the ITM for /3 = 



— 0. 15 and reverse flow iterations. 



/ 

J 




~fnil Pn 


C vlALl/HAlliJUO 





447 


47 


942 


1 
1 


^^^^ 


8/1 


1017 
IZl / 


2 


458 


47 


964 


3 


483 


47 


1014 


4 


464 


52 


981 


5 


425 


58 


909 


6 


479 


49 


1008 


7 


463 


46 


973 



Table 4: Comparison for the reverse flow skin-friction coefficients — t(O). For 
all cases we used /?o = 15 and h\ = 25. The iterations were, from top to bottom 
line: 8, 7, 9, 7, and 7. 



/3 


Stewartson [I49H 


Asaithambi [2J 


Auteri et al. [|6l 


ITM 


-0.025 


-0.074 






-0.074366 


-0.05 


-0.108 






-0.108271 


-0.1 


-0.141 


-0.140546 


-0.140546 


-0.140546 


-0.15 


-0.132 


-0.133421 


-0.133421 


-0.133421 


-0.18 


-0.097 


-0.097692 


-0.097692 


-0.097692 
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